Pseudo-potential treatment of two aligned dipoles under external harmonic confinement 
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Dipolar Bose and Fermi gases, which are currently being studied extensively experimentally and theoretically, 
interact through anisotropic, long-range potentials. Here, we replace the long-range potential by a zero-range 
pseudo-potential that simplifies the theoretical treatment of two dipolar particles in a harmonic trap. Our zero- 
range pseudo-potential description reproduces the energy spectrum of two dipoles interacting through a shape- 
dependent potential under external confinement very well, provided that sufficiently many partial waves are 
included, and readily leads to a classification scheme of the energy spectrum in terms of approximate angular 
momentum quantum numbers. The results may be directly relevant to the physics of dipolar gases loaded into 
optical lattices. 
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I. INTRODUCTION 

Many -body systems with dipolar interactions have attracted 
a lot of attention recently. Unlike the properties of ultra- 
cold atomic alkali vapors, which can be described to a very 
good approximation by a single scattering quantity (the s- 
wave scattering length), those of dipolar gases additionally 
depend on the dipole moment. This dipole moment can be 
magnetic, as in the case of atomic Cr [1, 2], or electric, as 
in the case of heteronuclear molecules such as OH Jj, 
KRb |5[] or RbCs |g]. Furthermore, dipolar interactions are 
long-ranged and anisotropic, giving rise to a host of novel 
many-body effects in confined dipolar gases such as roton- 



like features 



9 1 and rich stability diagrams ifiol [HI [H 
The physics of dipolar gases loaded into 
optical lattices promises to be particularly rich. For example, 
this setup constitutes the sta rting po int for a range of quan- 
tum computing schemes 11181 1 1 9L I20L Elf . Additionally, a va- 
riety of novel quantum phases have already been predicted to 
arise J22L I23L I24I EBIl . Currently, a number of experimental 
groups are working towards loading dipolar gases into optical 
lattices. 

This paper investigates the physics of doubly-occupied op- 
tical lattice sites in the regime where the tunneling between 
neighboring sites and the interactions with dipoles located in 
other lattice sites can be neglected. In this case, the prob- 
lem reduces to treating the interactions between two dipoles 
in a single lattice site. Assuming that the lattice potential 
can be approximated by a harmonic potential, the center of 
mass motion separates and the problem reduces to solving 
the Schrodinger equation for the relative distance vector r 
between the two dipoles. The interaction between the two 
aligned dipoles is angle-dependent and falls off as 1 /r 3 at 
large interparticle distances. In this work, we replace the 
shape-dependent interaction potential by an angle-dependent 
zero-range pseudo-potential, which is designed to reproduce 
the scattering properties of the full shape-dependent interac- 
tion potential, and derive an implicit eigenequation for two 
interacting identical bosonic dipoles and two interacting iden- 
tical fermionic dipoles analytically. 

Replacing the full interaction potential or a shape- 



dependent pseudo-potential by a zero-range pseudo- 
potential il HI El |H [M [M [H often allows for an 
analytical description of ultracold two-body systems in 
terms of a few key physical quantities. Here we show 
that the eigenequation for appropriately chosen zero-range 
pseudo-potentials reproduces the energy spectrum of two 
dipoles under harmonic confinement interacting through 
a shape-dependent model potential; that the applied zero- 
range treatment readily leads to an approximate classification 
scheme of the energy spectrum in terms of angular momentum 
quantum numbers; and that the proposed pseudo-potential 
treatment breaks down when the characteristic length of the 
dipolar interaction becomes comparable to the characteristic 
length of the external confinement. The detailed understand- 
ing of two interacting dipoles obtained in this paper will 
guide optical lattice experiments and the search for novel 
many-body effects. 

Section [TT] introduces the Hamiltonian under study and dis- 
cusses the anisotropic zero-range pseudo-potential that is used 
to describe the scattering between two interacting dipoles. In 
Sec. [in] we derive an implicit eigen equation for two dipoles 
under external spherical harmonic confinement interacting 
through the zero-range pseudo-potential and show that the 
resulting eigenenergies agree well with those obtained for a 
shape-dependent model potential. Finally, Sec. HVIconcludes. 



II. SYSTEM UNDER STUDY AND ANISOTROPIC 
PSEUDO-POTENTIAL 

Within the mean-field Gross-Pitaevskii formalism, the in- 
teraction between two identical bosonic dipoles, aligned along 
the space-fixed z-axis by an external field, has been success- 
fully modeled by the pseudo-potential V pp (r) ifllll . 
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(1) 



Here, /j denotes the reduced mass of the two-dipole system, 
d the dipole moment, and the angle between z and the rel- 
ative distance vector r. The s-wave scattering length aqo de- 



2 



pends on both the short- and long-range parts of the true in- 
teraction potential. The second term on the right hand side of 
Eq. ([]]) couples angular momentum states with I = I' (I > 0) 
and \l — l'\=2 (any /,/')■ F° r identical fermions, s-wave scat- 
tering is absent and the interaction is described, assuming the 
long-range dipole-dipole interaction is dominant, by the sec- 
ond term on the right hand side of Eq. ([1). 

Our goal in this paper is to determine the eigenequation 
of two identical bosonic dipoles and two identical fermionic 
dipoles under external spherically harmonic confinement with 
angular trapping frequency CO analytically. The Schrodinger 
equation for the relative position vector r reads 



with 



[H Q + V mt {r)]W)=EW), 



(2) 



where the Hamiltonian Ho of the non-interacting harmonic os- 
cillator is given by 



Ho 



2/j 



1 



2 2 



-/j to r 



(3) 



In Eq. (O, Vi„ t (r) denotes the interaction potential. The 
pseudo-potential V pp (r) cannot be used directly in Eq. (f2]i 
since both parts of the pseudo-potential lead to divergencies. 
The divergence of the 8-function potential arises from the sin- 
gular 1 jr behavior at small r of the spherical Neumann func- 
tion no(r), and can be cured by introducing the regularization 
operator l27ll . Curing the divergence of the long -ranged 
1 /r 3 term of V pp is more involved, since it couples an infi- 
nite number of angular momentum states, each of which gives 
rise to a singularity in the r — > limit. The nature of each of 
these singularities depends on the quantum numbers I and /' 
coupled by the pseudo-potential, and hence has to be cured 
separately for each / and /' combination. 

In this work, we follow Derevianko ll33il34ll and cure the di- 
vergencies by replacing V pp (r) with a regularized zero-range 
potential Vp p reg (r), which contains infinitely many terms, 



V n 



IV 



(4) 



The sum in Eq. © runs over I and /' even for identical bosons, 
and over I and I' odd for identical fermions. For I ^ I', Vu> 
and Vpi are different and both terms have to be included in 
the sum. In Sec. [HI] we apply the pseudo-potential to sys- 
tems under spherically symmetric external confinement. For 
these systems, the projection quantum number m is a good 
quantum number, i.e., the energy spectrum for two interact- 
ing dipoles under spherically symmetric confinement can be 
solved separately for each allowed to value. Consequently, a 
separate pseudo-potential can be constructed for each m value. 
In the following, we restrict ourselves to systems with vanish- 
ing projection quantum number to; the generalization of the 
pseudo-potential to general to is discussed at the end of this 
section. The V/// are defined through their action on an arbi- 
trary r-dependent function < t ) (r) 13 3l 13411 . 

V ll ,(r)<P(r)=g ll ,^Y, /0 (Q^)x 



d 2l+l 



3 r 2/+l 



J Y l0 (B,$)®(r)dQ, 



(5) 



gw 



h 2 a w (2Z+l)!!(2/' + l)! 



(2/ + 1)! 



(6) 



where k denotes the relative wave vector, k — y/ 2/jE / h 2 , and 
the aui generalized scattering lengths. Since we are restricting 
ourselves to to = 0, the V/// are written in terms of the spherical 
harmonics F/„, with to = 0. When applying the above pseudo- 
potential we treat a large number of terms in Eq. @, and do 
not terminate the sum after the first three terms as done in 
Refs. JUHB. We note that the non-Hermiticity of V P p^eg 
does not lead to problems when determining the energy spec- 
trum; however, great care has to be taken when calculating, 
e.g., structural expectation values 13611 . 

To understand the functional form of the zero-range 
pseudo-potential defined in Eqs. through ©, let us first 
consider the piece of Eq. Q in square brackets. If we decom- 
pose the incoming wave <I>(r) into partial waves, 



*(?)= L c n . limi Q nili (r)Y lm (e,ty), 



(7) 



njljnij 



where the c n; /. m . denote expansion coefficients and the Q ni i 
radial basis functions, the spherical harmonic F/o in the inte- 
grand of V/// acts as a projector or filter. After the integra- 
tion over the angles, only those components of < t ) (r) that have 



lj = I and to,- = survive. The operator J-j^ r /+1 in Eq. <(5j is 

designed to then first cure the r~'~' divergencies of the Q nj i, 
which arise in the r — » limit, and to then second "extract" the 
coefficients of the regular part of the Q nj i (r) that go as r l l27ll . 
Alltogether, this shows that the square bracket in Eq. © re- 
duces to a constant when the r — > limit is taken. To under- 
stand the remaining pieces of the pseudo-potential, we mul- 
tiply Eq. d5]l from the left with Q* , F,* and integrate over 

A * x — n o L o h m o j . ^ 

all space. The spherical harmonic F//o in Eq. (|5jl then ensures 
that the integral is only non-zero when /' = l (> and m () = 0. 
When performing the radial integration, the 8(r)/r' term en- 
sures that the coefficients of the regular part of the Q n „i that 
go as r l ° are being extracted (note that the remaining 1 jr 2 term 
cancels the r 2 in the volume element). 

Alltogether, the analysis outlined in the previous paragraph 
shows that the functional form of V\f ensures that the diver- 
gencies of the radial parts of the incoming and outgoing wave 
is cured in the r — * limit and that the /th component of the 
incoming wave is scattered into the /'th partial wave. The 
sum over all I and I' values in Eq. (0]l guarantees that any 
state with quantum number I can be coupled to any state with 
quantum number /', provided the corresponding generalized 
scattering length is non-zero. We note that the regularized 
pseudo-potential given by Eqs. © through © is only appro- 
priate if the external confining potential in Eq. (f3]) has spher- 
ical symmetry i37tl . Generalizations of the above zero-range 
pseudo-potential, aimed at treating interacting dipoles under 
elongated confinement, require the regularization scheme to 
be modified to additionally cure divergencies of cylindrically 
symmetric wave functions. These extensions will be subject 
of future studies. 
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We now discuss the generalized scattering lengths a///, 
which determine the scattering strengths of the V///. The a n < 
have units of length and are defined through the K-matrix ele- 

' 3, 



ments K\™ 



«//' 



lim ■ 



(8) 



for to — 0. The scattering lengths a/// and am are identical 
because the K-matrix is symmetric. In general, the scattering 
lengths a hi have to be determined from the K-matrix elements 
for the "true" interaction potential, which contains the long- 
range dipolar and a short-ranged repulsive part, of two inter- 
acting dipoles. As discussed further in Sec. [TTT1 an approach 
along these lines is used to obtain the squares shown in Fig. [3] 
Alternatively, it has been shown that the K-matrix elements 
(except for Kqq , see below) for realistic potentials, such as for 
the Rb-Rb potential in a strong electric field [11] or an OH- 
OH model potential [16], are approximated with high accu- 
racy by the K-matrix elements for the dipolar potential only, 
calculated in the first Born approximation. Applying the Born 
approximation to the second term on the right hand side of 
Eq. ©, we find for to = and I = /'(/> 1) 



2ZX 



an 



(21- l)(2Z + 3) : 
and for m = and I = I' + 2 



(2/-l) v /(2/ + l)(2/-3)' 



(9) 



(10) 



For I 1 = 2 and I — 0, e.g., Eq. ( fTOb reduces to 020 = 
— D*/(3\/5), in agreement with Ref. 113 311 . The scattering 
lengths fl/_2,/ are equal to a/ /_2> an d a U other generalized scat- 
tering lengths are zero. In Eqs. (O and ( fTOb , denotes the 
dipole length, = /jd 2 /h 2 . All non-zero scattering lengths 
aw are negative, depend on I and and are directly propor- 
tional to d 2 . Furthermore, for fixed D„, the absolute value of 
the non-zero a^ decreases with increasing angular momen- 
tum quantum number /, indicating that the coupling between 
different angular momentum channels decreases with increas- 
ing /. However, this decrease is quite slow and, in general, an 
accurate description of the two-dipole system requires that the 
convergence with increasing be assessed carefully. 

One can now show readily that the K-matrix elements K 1 ^ 
of V PPt regi calculated in the first Born approximation, with a^ 
given by Eqs. (O and ( fTOb coincide with the K-matrix ele- 
ments K 1 ^ of V pp . This provides a simple check of the zero- 
range pseudo-potential construction and proofs that the pref- 
actors of Vw are correct. In turn, this suggests that the appli- 
cability regimes of V pp and V PPl reg are comparable, if the gen- 
eralized scattering lengths a/// used to quantify the scattering 
strengths of Viji are approximated by Eqs. (O and ( fTOb . The 
applicability regime of V pp>reg may, however, be larger than 
that of V pp if the full energy-dependent K-matrix of a realistic 
potential is used instead. 

To generalize the zero-range pseudo-potential defined in 
Eqs. (O through (O for projection quantum numbers m = 



to any to, only a few changes have to be made. In Eq. ©, 
the spherical harmonics F/o have to be replaced by F/ m , and 
the generalized scattering lengths have to be defined through 
limfc_>Q — Kl™ jk. Correspondingly, Eqs. © and (TTOb become 
TO-dependent. 



III. TWO DIPOLES UNDER EXTERNAL CONFINEMENT 



Section IIII Al derives the implicit eigenequation for two 
dipoles interacting through the pseudo-potential under exter- 
nal harmonic confinement and Section IIII B I analyzes the re- 
sulting eigen spectrum. 



A. Derivation of the eigenequation 

To determine the eigen energies of two aligned dipoles 
with to = under spherical harmonic confinement interacting 
through the zero-range potential Vpp je g, we expand the eigen- 
functions *P(r) in terms of the orthonormal harmonic oscilla- 
tor eigen functions F/,o, 



¥(r)= 2>„ iZi /? nj/i .(r)F /i0 (e, 



(11) 



The pseudo-potential V PPire g enforces the proper boundary 
condition of *P(r) at r — 0, and thus determines the expan- 
sion coefficients c n; / ; . To introduce the key ideas we first con- 
sider i-wave interacting particles rf28tl . for which the pseudo- 
potential reduces to a single term, and then consider the gen- 
eral case, in which the pseudo-potential contains infinitely 
many terms. 

Including only the term with / and /' = in Eq. ©, the 
Schrodinger equation becomes, 

Y,c nik {E nil .-E + V m )R nik {r)Y li0 {Q,if) - 0, (12) 

no- 
where the E n .i denote the eigenenergies of the non-interacting 
harmonic oscillator, 



(13) 



E„ ih = [ 2m+lt + - ) to. 



In what follows, it is convenient to express the energy E of the 
interacting system in terms of a non-integer quantum number 
V, 



E = ( 2v + - ) ftco. 



(14) 



Multiplying Eq. ([T2]i from the left with R* n , Yf with l > 
and integrating over all space, we find that the c„ ( ./ ; with 
/, > vanish. This can be understood readily by realizing 
that the s-wave pseudo-potential Voo, as discussed in detail in 
Sec. HI] only couples states with I — I' = 0. To determine the 
expansion coefficients c n ,o, we multiply Eq. (fT2l from the left 
with ^* o oF *o and integrate over all space. This results in 

c no0 (2n o - 2v)m + R* no0 (Q)gooB - 0, (15) 
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where Bq denotes the result of the square bracket in Eq. (0, 



( r £ c ni0 R„fi(r) 



(16) 



Note that Bq is constant and independent of n,. In Eq. ( fT3T l, 
the r-independent term ^* o q(0) arises from the radial integra- 
tion over the 8-function of the pseudo-potential. If we solve 
Eq. ([TBI l for c, Io o and plug the result into Eq. ( TToT l. the unknown 
constant Bo cancels and we obtain an implicit eigenequation 
for v, 



goo 




^o(0)^o(0 ' 
n ^ (2v-2m)tm 



(17) 



Using Eqs. dZTb and d32b from the Appendix to simplify the 
term in square brackets, we obtain the well-known implicit 
eigenequation for two particles interacting through the s-wave 
pseudo-potential under spherical harmonic confinement [28], 



1 \2fUO 



2r 



' -E 
.. 2ho} 



— — o 

aim 



(18) 



H ere, ai w denotes the harmonic oscillator length, = 

-v/ft/Cu©). 

The derivation of the implicit eigenequation for two dipoles 
under external harmonic confinement interacting through the 
pseudo-potential with infinitely many terms proceeds analo- 
gously to that outlined above for the s-wave system. The key 
difference is that each V H / term in Eq. Q with I ^ I' couples 
states with different angular momenta, resulting in a set of 
coupled equations for the expansion coefficients c nj ; ; . How- 
ever, since V pp ^ reg for dipolar systems couples only angular 
momentum states with |/ — l'\ < 2 [see, e.g., the discussion at 
the beginning of Sec.|II]and around Eqs. (0 and (TTOlil. the cou- 
pled equations can, as we outline in the following, be solved 
analytically by including successively more terms in V PP}reg . 

To start with, we plug the expansion given in Eq. (fTTT i into 
Eq. (f2]l, where the interaction potential V inf is now taken to 
be the pseudo-potential V pp , reg with infinitely many terms. To 
obtain the general equation for the expansion coefficients c n; /j, 
we multiply as before from the left with R* , F,* and integrate 
over all space, 



c„ ij2n + l -2v)h(ti + 



\glo-2,l„Bl -2 + gl l Bl +gl B +2J B Bl 0+ 2] = 0. 



(19) 



Here, the B/ -2> B\ o and ^/„+2 denote constants that are inde- 
pendent of «,-, 



Bl a = 



^2l„ + l 



Jo+l 



E C n,l Rn,lo( r ) 
\"i=° J 



(20) 



J r^O 



The three terms in the square bracket in the second line of 
Eq. ( fT9l ) arise because the V/'_2/, Vw and Vl>+21> terms in the 
pseudo-potential V p 



'pp.reg couple the state R* t Y { * , for I' = l , 



with three components of the expansion for X V, Eq. ( TTTb . 
Importantly, the constants fi; -2, B[ o and B[ o+ 2, defined in 
Eq. ( f20b . depend on the quantum numbers l„ — 2, l and l + 2, 
respectively, which implies that Eq. ( fT9l ) defines a set of in- 
finitely many coupled equations that determine, together with 
Eq. ( 1201 ). the expansion coefficients c nj i r Notice that Eqs. ( fT9b 
and ( f20b coincide with Eqs. ( TT3T > and ( TToT ) if we set l = and 
= if I or /' > 0. 
We now illustrate how Eqs. ([T9i l and (f20b can be solved for 
identical bosons, i.e., in the case where I and Z' are even (the 
derivation for identical fermions proceeds analogously). Our 
strategy is to solve these equations by including successively 
more terms in the coupled equations, or equivalently, in the 
pseudo-potential. As discussed above, if aoo is the only non- 
zero scattering length, the eigenenergies are given by Eq. ( TT8l . 
Next, we also allow for non-zero «2o, «02 and 022, i.e., we con- 
sider I and /' < 2 in Eq. (@}. In this case, the coefficients c„ ; o 
and c nj 2 are non-zero and coupled, but all c„./. with /, > 2 are 
zero. Using the expressions for Bq and B2 given in Eq. d20l ), 
we decouple the equations. Finally, using Eqs. ( f27l > and (f32b 
from the Appendix, the eigenequation can be compactly writ- 
ten as 



where 



and 



r(— - 



to- 



■£=0, 



2) 



2 2Z+l r( '^£ . 
z 1 ( 2Bco 



«// 



k 2l a 2l+l 
K a ho 



<ll : 



2,/ 



. 4/-4„4i-2 ■ 
K "ho 



(21) 



(22) 



(23) 



Equation (l2Tb can be understood as follows. If only aoo is 
non-zero, it reduces to fo = 0, in agreement with Eq. ( fT8l . If 
only aoo, ao2 and 020 are non-zero, Eq. (ETT i remains valid if 
a22 in t2 is set to zero. This shows that the term q2 and the 
first term on the right hand side of t2 arise due to the coupling 
between states with angular momenta and 2. The second 
term of t-i, in contrast, arises due to a non-zero 022- Finally, 
for non-zero aoo and 022 but vanishing 020 and 002, Eq. ( l2Tb 
reduces to fo?2 = 0. In this case, we recover the eigenequations 
1 = for s-wave interacting particles j28tl and ti = for d- 
wave interacting particles IB2I1 . 

We now consider I and I' values with up to l max = 4 in 
Eq. (01, i.e., we additionally allow for non-zero 024, «42 and 
«44, and discuss how the solution changes compared to the 
lmax = 2 case. The equation for the expansion coefficients 
c m Q remains unchanged while that for c„,2 is modified. Fur- 
thermore, the expansion coefficients c Bj 4 are no longer zero. 
Consequently, we have three coupled equations, which can be 
decoupled, resulting in the following implicit eigenequation, 
to + ml (t2 + /h) = 0. In analogy to the l max = 2 case, the q\ 
term and the first part on the right hand side of the t\ term arise 
due to the "off-diagonal" scattering lengths 024 and 042, and 
the second term of t\ arises due to the "diagonal" scattering 
length 044. 
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Next, let us assume that we have found the implicit 
eigenequation for the case where we include terms in Eq. © 
with / and /' up to l max — 2. If we now include terms with 
/ and I 1 up to l max , only the equations for the expansion co- 
efficients c„ i B with /„ = l max — 2 and l max change; those for 
the expansion coefficients c„ / with l < l max — 4 remain un- 
changed. This allows the l max /2 + 1 coupled equations for the 
expansion coefficients to be decoupled analytically using the 
results already determined for the case where / and I' go up 
to lmax — 2. Following this procedure, we find the following 
implicit eigenequation 



T, = 0, 

'max ' 



(24) 



where 7/ mav itself can be written as a continued fraction. For 
identical bosons we find, 



to- 



qi 



<74 



(25) 



Taking l max — > °° gives the eigenequation for two identi- 
cal bosons under spherical harmonic confinement interacting 
through Vppj- e g with infinitely many terms. For two identical 
fermions, Eqs. ( 1221 through d25l l remain valid if the subscripts 
0,2, • • • in Eq. d25l l are replaced by 1 ,3, ■ • •. 

The derived eigenequation reproduces the eigenenergies in 
the known limits. For the non-interacting case (all a/// = 0), 
the eigenenergies coincide with the eigenenergies of the har- 
monic oscillator, i.e., E„i — (2n + I + 3/2)fixo, where n = 
0, 1,2, ■ • • and / — 0,2,4, ■ ■ ■ (in the case of identical bosons) 
and I = 1 , 3 , ■ • • (in the case of identical fermions). The kth lev- 
els, with energy (2k + 3/2)h(ti for bosons and (2fe + 5/2)fi,co 
for fermions, has a degeneracy of k+ 1, k = 0, 1, • • •. Non- 
vanishing a/;/ lead to a splitting of degenerate energy levels 
but leave the number of energy levels unchanged. If a// is 
the only non-zero scattering length, the eigenequation reduces 
to that obtained for spherically symmetric pseudo-potentials 
with partial wave / 03211 . 



B. Analysis of the energy spectrum 

This section analyses the implicit eigenequation, Eq. (1241 . 
derived in the previous section for the zero-range pseudo- 
potential for m = and compares the resulting energy spec- 
trum with that obtained for a shape-dependent model poten- 
tial. The implicit eigenequation, Eq. ( |24"1 >. can be solved 
readily numerically by finding its roots in different energy 
regions. The solutions of the Schrodinger equation for the 
shape-dependent model potential are otained by expanding the 
eigenfunctions on a B-spline basis. 

Lines in Figs.QJa) and (b) show the eigenenergies obtained 
by solving Eq. (1241 for two identical bosons and two iden- 
tical fermions, respectively, interacting through Vp S , re g(r) un- 
der external spherically symmetric harmonic confinement as 
a function of the dipole length D*. In both panels, we as- 
sume that the interaction between the two dipoles is purely 
dipolar, i.e., in Fig.[Tfa) we set aoo = 0. The other scattering 
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Figure 1: Relative eigenenergies E for (a) two identical bosonic 
dipoles and (b) two identical fermionic dipoles interacting through 
Vpp,reg [using oqo = in (a)] under spherical harmonic confinement 
as a function of £>*/a/, . The line style indicates the predominant 
character of the corresponding eigenstates. In (a), a solid line refers 
to I ~ 0, a dashed line to I ~ 2, and a dotted line to / ~ 4; in (b), a 
solid line refers to / ~ 1, a dashed line to I ~ 3, and a dotted line to 
1^5. 



lengths a/// are approximated by Eqs. (O and ( fTOb . Interest- 
ingly, for identical bosons, the lowest gas-like level, which 
starts at E — 1 .5S.C0 for D„ = 0, increases with increasing Z)„. 
For identical fermions, in contrast, the lowest gas-like state 
decreases with increasing . 

In addition to obtaining the eigenenergies themselves, the 
pseudo-potential treatment allows the spectrum to be classi- 
fied in terms of angular momentum quantum numbers. To 
this end, we solve the implicit eigenequation, Eq. ( |24| ). for 
increasing l max , and monitor how the energy levels shift as ad- 
ditional angular momenta are included in V P p^ eg . Since a level 
with approximate quantum number / changes only little as 
larger angular momentum values are included in the pseudo- 
potential, this analysis reveals the predominant character of 
each energy level. In Fig. [If a), the eigenfunctions of energies 
shown by solid, dashed and dotted lines have predominantly 
/ = 0, 2 and 4 character, respectively. In Fig.[Tfb), the eigen- 
functions of energies shown by solid, dashed and dotted lines 
have predominantly / = 1 , 3 and 5 character, respectively. We 
find that the lowest excitation frequency between states with 
predominantly I = [/ = 1] character, increases [decreases] 
for identical bosons [fermions] with increasing . These pre- 
dictions can be verified directly experimentally. 

To assess the accuracy of the developed zero-range pseudo- 
potential treatment, we consider two interacting bosons with 
non-vanishing s-wave scattering length aoo- We imagine that 
the dipole moment of two identical polarized bosonic polar 
molecules is tuned by an external electric field. As the dipole 
moment d is tuned, the s-wave scattering length aoo, which 
depends on the short-range and the long-range physics of the 
"true" interaction potential, changes. To model this situation, 
we solve the two-body Schrodinger equation, Eq. (0, numer- 
ically for a shape-dependent model potential with hardcore 
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Figure 2: Panel (a) shows the relative energies E for two aligned 
identical bosonic dipoles under external spherical harmonic confine- 
ment as a function of £)* / a/ w . Solid lines show the numerically deter- 
mined energies obtained using V„ w del with b = 0.0097a/ lo . Crosses 
show the energies obtained using Vpp^reg with essentially infinitely 
many terms, and aoo calculated for V mot i e i- Panel (b) shows a blow- 
up of the energy region around E ~ 5.5fico. Note that the horizontal 
axis in (a) and (b) are identical. 



radius b and long-range dipolar tail. In this case, V,„ r is given 

by 



f ^i^We ifr>b 

Vmodelir) = < ? t Z U 

\ oo if r < b 



(26) 



For d = 0, the s-wave scattering length aoo for V mo dei is given 
by b. As the dipole length Z)„ increases, aoo g° es through 
zero, and becomes negative. Just when the two-body poten- 
tial supports a new bound state, aoo goes through a resonance 
and becomes large and positive. As increases further, aoo 
decreases. This resonance structure repeats itself with increas- 
ing D* (see Fig. 1 of Ref. [17]; note, however, that the lengths 
fl/io and defined throughout the present work differ from 
those defined in Ref. lUzl]). 

For the model potential V„ w deh aoo depends on the ratio 
between the short-range and long-range length scales, i.e., 
on b/D*. To compare the pseudo-potential energies and the 
energies for the model potential, we fix b and calculate aoo 
for each D t . considered. The dipole-dependent s-wave scat- 
tering length is then used in the zero-range pseudo-potential 
Vpp,reg- The other scattering lengths are, as before, approx- 
imated by the expressions given in Eqs. (0 and (fTOb - Solid 
lines in Fig. [2ja) and (b) show the eigenenergies obtained for 
Vmodei as a function of . Crosses show the eigenenergies ob- 
tained for Vpp. re g using a value of l max that results in converged 
eigenenergies. The overview spectrum shown in Fig. |2{a) 
shows that one of the energy levels dives down to negative 
energies close to that value at which the two-body poten- 
tial V nw dei supports a new bound state. The blow-up, Fig.|2jb), 
around E w 5.5?iCO shows excellent agreement between the en- 
ergies obtained using V P p. reg (crosses) and those obtained us- 
ing V mo dei (solid lines); the maximum deviation for the energy 



range shown is 0.05 %. 

As before, we can assign approximate quantum numbers to 
each energy level. At D* <C a/, , the three energy levels around 
E w 5.5h(£> have, from bottom to top, approximate quantum 
numbers I =2, A and 0. After two closely spaced avoided 
crossings around w 0.025a/,„, the assignment changes to 
I = 0, 2 and 4 (again, from bottom to top). If the maximum 
angular momentum l max of the pseudo-potential is set to 2, the 
energy level with approximate quantum number I — 4 would 
be absent entirely. This illustrates that a complete and accu- 
rate description of the energy spectrum requires the use of a 
zero-range pseudo-potential with infinitely many terms. The 
energy of a state with approximate quantum number / requires 
lmax to be at least I for the correct degeneracy be obtained and 
at least I + 2 for a quantitative description. 

The sequence of avoided crossings at « 0.025a/ !O sug- 
gests an interesting experiment. Assume that the system 
is initially, at small electric field (i.e., small /«/,„), pre- 
pared in the excited state with angular momentum I i=s and 
E ps 5.52^co. The electric field is then slowly swept across the 
first broad avoided crossing at 0.019a/ !O to transfer the 
population from the state with I « to the state with / « 2. 
We then suggest to sweep quickly across the second narrower 
avoided crossing at D* w 0.028a/, o (the ramp speed must be 
chosen so minimize population transfer from the state with 
/ « 2 to the state with I ss 4). As in the case of s-wave scatter- 
ing only (HI, the time-dependent field sequence has to be op- 
timized to obtain maximal population transfer. The proposed 
scheme promises to provide an efficient means for the transfer 
of population between states with different angular momenta 
and for quantum state engineering. 

Figure [2] illustrates that the pseudo-potential treatment re- 
produces the eigenenergies of the shape-dependent model po- 
tential Vmodei- To further assess the validity of the pseudo- 
potential treatment, we now consider two interacting bosonic 
dipoles for which the dipolar interaction is dominant, i.e., we 
consider aoo = 0. For V mo dei with b = 0.003 la/, , we deter- 
mine a set of values at which aoo = 0. Note that the num- 
ber of bound states with predominantly s-wave character in- 
creases by one for each successively larger D*. Crosses in 
Figs. HJa)-(c) show the eigenenergies for V mo dei with aoo = 
as a function of in the energy ranges around 1.5, 3.5 and 
5.5h(x). For comparison, lines show the eigenenergies ob- 
tained for the regularized pseudo-potential with aoo = 0. As in 
Fig.[T] the linestyle indicates the predominant character of the 
energy levels (solid line: I « 0; dashed line: / « 2; and dotted 
line: I « 4). The agreement between the energies obtained for 
the pseudo-potential with a hi given by Eqs. (0 and ( fTOb and 
for the model potential for small D„ is very good, thus vali- 
dating the applicability of the pseudo-potential treatment. The 
agreement becomes less good, however, as D„ increases. This 
can be explained readily by realizing that the dipole length Z)„ 
approaches the harmonic oscillator length a\ w - 

In general, the description of confined particles interacting 
through zero-range pseudo-potentials is justified if the char- 
acteristic lengths of the two-body potential are smaller than 
the characteristic length of the confining potential. For exam- 
ple, in the case of s-wave interactions only, the van der Waals 
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Figure 3: Crosses show the relative eigenenergies £ as a function 
of D*/a/ m for two identical bosons with aoo = interacting through 
Vmodel with b = 0.003 la^g in three different energy regions. Lines 
show E for two identical bosons with aoo = interacting through 
Vps.reg with anf given by Eq. ((9} and JlOb . As in Fig. fjja) solid, 
dashed and dotted lines show the energies of levels characterized by 
approximate quantum numbers / ~ 0, 2 and 4. The agreement be- 
tween the crosses and the lines is good at small £>„/a/, but less good 
at larger O*/a/, . Squares show the eigenenergies obtained for the 
energy-dependent pseudo-potential at £>» = 0.242a/,,,; the agreement 
between the squares and the crosses is excellent, illustrating that us- 
age of the energy-dependent K-matrix greatly enhances the applica- 
bility regime of V pp , re g ■ 

length has to be smaller than the oscillator length j29l l30ll . 
The model potential V mot i e i is characterized by a short-range 
length scale, the hardcore radius b, and the dipole length 
D*; in Fig. [3] it is the relatively large value of D*/ai w that 
leads, eventually, to a break-down of the pseudo-potential 
treatment. As in the case of spherical interactions, the break- 
down can be pushed to larger values by introducing 
energy-dependent generalized scattering lengths de- 
fined through —KiqQc) jk for m — 0, and by then solving the 
eigenequation, Eq. ( E4b . self-consistently J29l[30ll . 

Figure [4] shows three selected scattering lengths ani(k) for 
the model potential V 'model with = 18.9b as a function of 
energy. This two-body potential supports eight bound states 
with projection quantum number m = 0, which have predom- 
inantly s-wave character. Both energy and length in Fig.|4]are 
expressed in oscillator units to allow for direct comparison 
with the data shown in Fig. [3] The scattering length aoo(^), 
shown by a solid line in Fig. [4] is zero at zero energy and in- 
creases with increasing energy. Both a2o(k) (dashed line) and 
«22 (k) (dash-dotted line) are negative. Their zero-energy val- 
ues coincide with those calculated in the Born approximation 
(horizontal dotted lines). 

Using these energy-dependent «//'(£) to parametrize the 
strengths of the pseudo-potential and solving the eigenequa- 
tion, Eq. d24l ). self-consistently, we obtain the squares in 
Fig. [3] The energies for V pp ^ e g with energy-dependent a[[i 
(squares) are in much better agreement with the energies ob- 
tained for the model potential (crosses) than the energies ob- 
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Figure 4: Energy-dependent scattering lengths am(k) (solid line), 
a 2o(k) (dashed line) and axi{k) (dash-dotted line) for the model po- 
tential V moc iei with D t = 78.9fc as a function of the relative energy 
E. In oscillator units, V mo j e i is characterized by b = 0.003 la/ 1( , and 
£>* = 0.242a/ lo . For comparison, horizontal dotted lines show the 
energy-independent scattering lengths «22> Eq. (9^, and 020, Eq. dlOt . 
calculated in the first Born approximation. 



tained using the energy-independent ajji to parametrize the 
pseudo-potential (lines). This suggests that the applicability 
regime of the regularized zero-range pseudo-potential can be 
extended significantly by introducing energy-dependent scat- 
tering lengths. Since the proper treatment of resonant inter- 
actions within the regularized zero-range pseudo-potential re- 
quires that the energy-dependence of the generalized scatter- 
ing lengths be included, future work will address this issue in 
more depth. 



IV. SUMMARY 

This paper applies a zero-range pseudo-potential treatment 
to describe two interacting dipoles under external spherically 
harmonic confinement. Section [XT] introduces the regular- 
ized zero-range pseudo-potential V pp ^ reg used in this work, 
which was first proposed by Derevianko 133l 13411 . Particu- 
lar emphasis is put on developing a simple interpretation of 
the individual pieces of the pseudo-potential. Furthermore, 
we clearly establish the connection between V ppje g and the 
pseudo-potential V pp , which is typically employed within a 
mean-field framework. We argue that the applicability regime 
of these two pseudo-potentials is comparable if the scattering 
strengths of V pPtreg , calculated in the first Born approximation, 
are chosen so as to reproduce those of V pp . 

We then use the regularized zero-range pseudo-potential to 
derive an implicit eigen equation for two dipoles under ex- 
ternal confinement, a system which can be realized experi- 
mentally with the aid of optical lattices. In deriving the im- 
plicit eigenequation, we again put emphasis on a detailed 
understanding of how the solution arises, thus developing a 
greater understanding of the underlying physics. The im- 
plicit eigenequation can be solved straightforwardly, and al- 
lows for a direct classification scheme of the resulting eigen- 
spectrum. By additionally calculating the eigen energies for 
two dipoles interacting through a finite range model poten- 
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tial numerically, we assess the applicability of the developed 
zero-range pseudo-potential treatment. We find good agree- 
ment between the two sets of eigenenergies for small D* and 
quantify the deviations as increases. Finally, we show that 
the validity regime of V pp , reg can be extended by parametriz- 
ing the scattering strengths of V PPt reg in terms of the energy- 
dependent K-matrix calculated for a realistic model potential. 
This may prove useful also when describing resonantly inter- 
acting dipoles. 

At first sight it may seem counterintuitive to replace 
the long-range dipolar interaction by a zero-range pseudo- 
potential. However, if the length scales of the interaction po- 
tential, i.e., the van der Waals length scale characterizing the 
short-range part and the dipole length characterizing the long- 
range part of the potential, are smaller than the characteristic 
length of the trap ai w , this approach is justified since the zero- 
range pseudo-potential is designed to reproduce the K-matrix 
elements of the "true" interaction potential. This is particu- 
larly true if the pseudo-potential is taken to contain infinitely 
many terms, as done in this work. 

In summary, this paper determines the eigenenergies of two 
interacting dipoles with projection quantum number m = 0. 
The applied zero-range pseudo-potential treatment is vali- 
dated by comparing the resulting eigenenergies with those 
obtained numerically for a shape-dependent model potential. 
The analysis presented sheds further light on the intricate 
properties of angle-dependent scattering processes and their 
description through a regularized zero-range pseudo-potential 
with infinitely many terms. The calculated energy spectrum 
may aid on-going experiments on dipolar Bose and Fermi 
gases. 

Acknowledgements: KK and DB acknowledge support by 
the NSF through grant PHY-0555316 and JLB by the DOE. 



V. APPENDIX 



Writing the radial harmonic oscillator functions R n i(r) in 



terms of the Laguerre polynomials L„ 

Rnl{r) = 



(Z+l/2) 



2 l+2 



(2l + l)Un^ +1/2 \o)c 



ho 



-i-4)fe) >/2,(r2/ " L) ' ,2s) 



we find 



Rnl{r) 



2 , +2L (/+l/2) (0) 

Using Eqs. d28l ) and ( 129) , the C; can be rewritten as 

2 z+i 



(29) 



Q 



(2/ + l)!!7ii/2 fl 2'+ 3 



(/+1/2) ( (J _s2 



2ai 

ho 



to (v-n-|)to 



(30) 



We evaluate the infinite sum in Eq. ( f30b using the properties 
of the generating function [40], 



E 

n=0 



L^I 2 \{r/a ho f) 
v-n-i 



-r(-V + l/2)U(-V + 1/2,1 + 3/2, (r/a ho ) 2 ). (31) 

Using Eq. (l3~TT l together with the small r behavior of the hyper- 
geometric function U 114 111 , the expression for the Q reduces 
to 



In this Appendix, we evaluate the following infinite sum, 



Q 



j+i 



KM 



dr^ ^2{v-n-L) h& 



(27) 



Q 



(_l)/ 2 2/ + 2 (2/)!! £(-V+|) 1 
(2/ + 1)!! r(-v-^i)ftco^- 



(32) 
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